home *** CD-ROM | disk | FTP | other *** search
/ Collection of Tools & Utilities / Collection of Tools and Utilities.iso / ada / gnat1792.zip / gnat179b / t-adainc / s-xp_nbm.adb < prev    next >
Text File  |  1994-05-19  |  3KB  |  70 lines
  1. ------------------------------------------------------------------------------
  2. --                                                                          --
  3. --                         GNAT RUNTIME COMPONENTS                          --
  4. --                                                                          --
  5. --                        S Y S T E M . X P _ N B M                         --
  6. --                                                                          --
  7. --                                 B o d y                                  --
  8. --                                                                          --
  9. --                            $Revision: 1.2 $                              --
  10. --                                                                          --
  11. --           Copyright (c) 1992,1993,1994 NYU, All Rights Reserved          --
  12. --                                                                          --
  13. -- GNAT is free software;  you can  redistribute it  and/or modify it under --
  14. -- terms of the  GNU General Public License as published  by the Free Soft- --
  15. -- ware  Foundation;  either version 2,  or (at your option) any later ver- --
  16. -- sion.  GNAT is distributed in the hope that it will be useful, but WITH- --
  17. -- OUT ANY WARRANTY;  without even the  implied warranty of MERCHANTABILITY --
  18. -- or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License --
  19. -- for  more details.  You should have  received  a copy of the GNU General --
  20. -- Public License  distributed with GNAT;  see file COPYING.  If not, write --
  21. -- to the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. --
  22. --                                                                          --
  23. ------------------------------------------------------------------------------
  24.  
  25. function System.Xp_NBM
  26.   (Left    : Integer;
  27.    Modulus : Integer;
  28.    Right   : Natural)
  29.    return    Integer
  30. is
  31.    Result : Integer := 1;
  32.    Factor : Integer := Left;
  33.    Exp    : Natural := Right;
  34.  
  35.    function Mult (X, Y : Integer) return Integer;
  36.    pragma Inline (Mult);
  37.    --  Modular multiplication. Note that we can't take advantage of the
  38.    --  compiler's circuit, because the modulus is not known statically.
  39.  
  40.    function Mult (X, Y : Integer) return Integer is
  41.    begin
  42.       return Integer
  43.         (Long_Long_Integer (X) * Long_Long_Integer (Y)
  44.           mod Long_Long_Integer (Modulus));
  45.    end Mult;
  46.  
  47. --  Start of processing for System.Xp_NBM
  48.  
  49. begin
  50.    --  We use the standard logarithmic approach, Exp gets shifted right
  51.    --  testing successive low order bits and Factor is the value of the
  52.    --  base raised to the next power of 2.
  53.  
  54.    --  Note: it is not worth special casing the cases of base values -1,0,+1
  55.    --  since the expander does this when the base is a literal, and other
  56.    --  cases will be extremely rare.
  57.  
  58.    while Exp /= 0 loop
  59.       if Exp rem 2 /= 0 then
  60.          Result := Mult (Result, Factor);
  61.       end if;
  62.  
  63.       Factor := Mult (Factor, Factor);
  64.       Exp := Exp / 2;
  65.    end loop;
  66.  
  67.    return Result;
  68.  
  69. end System.Xp_NBM;
  70.